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Abstract 

We consider a particle system of the squared Bessel processes with index v > — 1 
conditioned never to collide with each other, in which if — 1 < v < the origin is 
assumed to be reflecting. When the number of particles is finite, we prove for any 
fixed initial configuration that this noncolliding diffusion process is determinantal in 
the sense that any multitime correlation function is given by a determinant with a con- 
tinuous kernel called the correlation kernel. When the number of particles is infinite, 
we give sufficient conditions for initial configurations so that the system is well de- 
fined. There the process with an infinite number of particles is determinantal and the 
correlation kernel is expressed using an entire function represented by the Weierstrass 
canonical product, whose zeros on the positive part of the real axis are given by the 
particle-positions in the initial configuration. Prom the class of infinite-particle initial 
configurations satisfying our conditions, we report one example in detail, which is a 
fixed configuration such that every point of the square of positive zero of the Bessel 
function J u is occupied by one particle. The process starting from this initial configu- 
ration shows a relaxation phenomenon converging to the stationary process, which is 
determinantal with the extended Bessel kernel, in the long-term limit. 
Keywords Noncolliding diffusion process, Squared Bessel process, Fredholm determi- 
nants, Entire functions, Weierstrass canonical products, Infinite particle systems 

1 Introduction 

Let 9JT be the space of nonnegative integer-valued Radon measures on R, which is a Polish 
space with the vague topology. We say £ n G DJl,n G N = {1,2,...} converges to £ G Wl 
weakly in the vague topology, if lim^oo f R ip{x)^ n {dx) = f R ip(x)£,(dx) for any ip G C (R), 
where Co(R) is the set of all continuous real- valued functions with compact supports in 
R. Any element of Wl can be represented by X^eA^O with a sequence of points in R, 
x = (xj) ie A) satisfying : x« G /} < oo for any compact subset / C R, and with a 
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countable set (an index set) A. We call an element £ of Wl an unlabeled configuration, 
and a sequence of points x a labeled configuration. For A C JR., we write the restriction 
of £ G 971 on A as (f n A)(-) = EieA^eA M - )- Let M + = {x G R : i > 0} and define 
9Jt + = {(£ nR+)(-) : £(•) G 371} . In the present paper we consider a one-parameter family of 
9JT + - valued processes with a parameter i/ > — 1, 

H M (*,0 = E^(*)(0. *e[o,oo), (i.i) 

describing a particle system of squared Bessel processes with index u > —1 (BESQ^) 
interacting with each other by long-ranged repulsive forces, such that X^(£)'s satisfy the 
SDEs 

dX\ u \t) = 2^jxl v \t)dB i {t) + 2{v + l)dt 

+ 4Xf } (t)y— r^dt, i = 1, 2, • ■ ■ , te[0,oo) (1.2) 

with a collection of independent standard Brownian motions (BMs), {Bi(t),i G N}, and, 
if — 1 < v < 0, with a reflection wall at the origin. Note that for the BM in R d , B(t) = 
(Bx(t), . . . , Bd(t)), d G N, the square of its distance from the origin, X(t) = \B(t)\ 2 = 
Eti Bi(t) 2 , solves the SDE, dX{t) = 2y/X(j)dB(t) + 2(v + l)dt with 

i/=^-l, (1-3) 

where -B(t) is a standard BM which is different from B^t), 1 < z < d [MJ H]. We give the 
initial configuration of the process £(•) = S^(0, •) = ^ &„.(•) and the process is denoted by 
(H<">(t),I*). 

When the number of particles is finite, £at(IR + ) = A^ < oo, the process (H^(t), P^) is 
realized in the following systems. 

(i) When v G No = N U {0}, i.e., when the corresponding dimension d given by (11.31) is a 
positive even integer, (S^(£), P^) is realized as the eigenvalue process of the Laguerre 
process [20]. Let M(t) be an (N + u) x N matrix, whose entries are independent 
complex BMs having the real and imaginary parts given by independent standard 
BMs, and set L(t) = M(t)*M(t). The N x N matrix-valued process L = (L(t)) te[0iOo) 
is called the Laguerre process. The matrix L(t) is Hermitian and positive definite, 
and its A" eigenvalues satisfy (ll.2p with % = 1, 2, . . . , N. When the entries of M(t) 
are independent standard real BMs, the matrix-valued process (M(t) T M(t)) tg [ 0oo ) is 
called the Wishart process [5], and thus L is also called the complex Wishart process. 
The eigenvalue processes of the real and complex Wishart processes are related with 
the random matrix theory [231 E] f° r the chiral Gaussian ensembles studied in the high 
energy physics (see [15] and references therein). 
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(ii) Let H(N) be the space of A^ x A" Hermitian matrices. And let sp (2 N, C) and so (2 N, C) 
be the symplectic Lie algebra and the orthogonal Lie algebra, having 2A" x 2A^-matrix 
representations, respectively. If the 2Nx2N matrix is in the space Hci^N) = H(2N)n 
sp(2N, C) or in U D (2N) = U(2N) Hso(2N, C), its eigenvalues are given by N pairs of 
positive and negative ones with the same absolute value, {(Aj, — Aj) : Aj > 0, 1 < i < 
N}. Consider the "Hc(2AQ-valued and the "HD(2AQ-valued Brownian motions. The 
dynamics of positive eigenvalues of them are described by (11.21) with v = 1/2 (d = 3) 
for the former case and with v = — 1/2 (d = 1) for the latter case, respectively [T5] . 
The pairing of positive and negative eigenvalues simulates the particle-hole symmetry 
in the energy space of the Bogoliubov-de Gennes formalism of superconductivity and 
these processes are related with the random matrix theory studied in the solid-state 
physics [Tj. 

(hi) Let p^(t, y\x), y e R+, be the transition probability density for BESQ^, v > —1, 



' i (y\ vj2 



P M (t,y\x) 



2t \x 



exp 



(2ty+ l T(v 
5(y - x), 



x + y 
2t 

~y/2t 



t>0,x>0, 
t > 0,x = 0, 



;i.4) 



if — 1 < v < 0, the origin is assumed to be reflecting 
Bessel function of the first kind defined by 



t = 0,x eR+, 
0], where I u (x) is the modified 



i»(x) = J2 



i 



n=0 



T(n+l)T(n+l+u) \2J 



x \ 2n + v 



1.5) 



with the Gamma function V{z) = / °° e u u z 1 du, > 0. First we consider the 
following Karlin-McGregor determinant 



f%\t,y\x)= det pMfayilxj) 

l<t,j<N L 



*>0, x = (x i )l 1 ,y = (y i )t 1 eW+, (1.6) 



where W 



N 



{X 



■ xn) e 



: < x 1 < 



< xjy}. The transition 



probability density of the A^-particle system of BESQ^ conditioned never to collide 
with each other, which we call the noncolliding BESQ^ U \ is given by the /i-transform 

of (USD, 

(1.7) 



pP(t,y\x) = h N (y)f^'(t,y\x 



h N (x) 



with the harmonic function given by the Vandermonde determinant [TTj (2D1 EES] 



h N (x) 



det \x] 

l<i,j<N 



n 



t 0C 4 



;i.8) 



l<i<j<N 
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It is easy to confirm that p^\t, -\x) satisfies the following backward Kolmogorov equa- 
tion 

d N d 2 - d 

g:u{t,x) = 2^Xi^u{t,x) + 2(v + l)^—u{t,x) 

i=l 1 i=l * 



N N r. 

Xi o 



i=l j=l J 



and it implies that the process (S^(t), P^), v > —1, is realized as the noncolliding 
BESQ M [E]. We put 

2JI+ = {£ g Wl + : £({x}) < 1 for any x G R}. (1.9) 

We see that SM(t) G 9J#, Vt > 0. 

In the present paper, we call the process (S ^(t),P$) the noncolliding BESQM. See 
, dni ET] for related noncolliding diffusion processes. 
Assume that £jy G VJt + with £jy(R + ) =iVeI For any M G N and any time sequence 
< t\ < ■ ■ ■ < tM < oo, the formula (11. 7p and the Markov property of the system give the 
multitime probability density of (H^P^) as [HI E] 



M ~ l f( u h+ (1)1 ^ 

^(ti,4 1} ;-..;tM,d M) ) = M* (M) ) J] f%\t m+1 -t m ;xW\xW y N {t h ^ x) |x) (1.10) 



m=l 



with £jv(") = Sj = i^i( - ) £ < xi < X2 < ■ ■ ■ < xn for the initial configuration and 

€n (") = Sili ^ .< m ) (') e a;(m ' ) = ( x i > • • • ? ^at ) e for the configurations at times 

tmi I < m < M. In fll.lOp . if some of x^s in x coincide, the factor ffc (ti, x^\x)/h N (x) is 
interpreted using l'Hopital's rule. 

For asM = {xt\ • • • , a#°) G W+ with = £^ * (m) (•) and N' G {1,2,..., iV}, 

we put cc]^ = (xi, . . . ,xffi) G W^,, 1 < m < M. For a sequence (iV m )^f =1 of positive 
integers less than or equal to N, we define the (N\, . . . , Nm)- multitime correlation function 
by 



'JVi' • • • ' fc Af i Nm 

M N M 



= L ,„._n n ^(*„^... i ««.d"')n w ^g f ,(i.n, 

,/ llm=i JK + m=li=AT m +l m=l v 7 

which is symmetric in the sense that pl N (- ■ ■ ; t m , a{xj^)\ . . . ) = p^ / iV (. . . ; t m , xffi] . . . ) with 

a(xj^) = (^^(l)' • • • ' ^(at™)) ^ or an ^ P ermu tation a G <Sjv m , 1 < Vm < M. For any M G N, 
/m G C (M+), m 6 M, 1 < m < M, < ti < • • • < t M < oo, the Laplace transform of (ll.lOp 
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is considered as a functional of x( x ) — {Xi{ x )i • • ■ j Xm( x ))-, where Xm( x ) = e e,n ^ m ^ — 1, 1 < 
m<M,xe K+; 



Ql N \x] = / n^ (ro) p^(ti,eSV--;^,d M) )ex P i^^ / / m (x)4 m) (^) 

^ K + M m=l * lrrt=l ^ M 

f M f 

U=i 



1.12) 



It is the generating function of multitime correlation functions, since if we expand it with 
respect to Xm(4™')' s ) pl N ' s appear as coefficients in terms; 



N M 



Gt N [x] 



AT , I ttX 7V m 



e "En? 

Af 1= N M =0m=l r 
M N m 

x II n^™^" 1 ^^ (*!' ^S; • • • ;*m, a; 



(M) 
% 



;i.i3) 



In the present paper, first we prove that, for any fixed initial configuration £jy £ 
with £tv(1R + ) = iV£N, there is a function IK-^s, x; i, ?/), which is continuous with respect 
to (x, y) G (0, oo) 2 for any fixed (s, t) G [0, oo) 2 , and that the function (I1.12p is given by the 
Fredholm determinant in the form 



u N lx\ 



Det 

l<m,n<M 



5 mn S{x -y) + Kl N (t m ,x;t n ,y)xn(y) 



;i.i4) 



By definition of Fredholm determinant (see Eq. fl4.8p in Section 4.1), (11.141) means that any 
multitime correlation function is given by a determinant 



(M)\ 
N M ) 



det 

l<i<N m ,l<j<N n 
Km,n<M 



mi •t'i 



;i.i5) 



The function "K\ N is called the correlation kernel and it determines the finite dimensional 
distributions of the process (E^(t),Fl N ) through (11.151) . It is an extension of determinantal 
(Fermion) point process of distributions studied by Soshnikov and Shirai and Takahashi 
[25] to the cases on T x R + with T = {t u . . . , t M }, M G N, < t x < ■ ■ ■ < t M < oo. See [12] 
for variety of examples of determinantal point processes. We express this result by simply 
saying that the noncolliding BESQ^ is determinantal with a correlation kernel K.l N for any 
6v e Wl + with &v( K +) = N G N (Theorem [2H]). 

Next we consider the infinite-particle limits. For £ G 9Jt + with = oo, when 

jg-?n[o,L] conver g eg ^ Q a continuous function as L — > oo, the limit is written as K£. If 
su Px, y ei |IK^ n [°' L l(s, x; t, y)| < oo,VL > for any (s,t) G (0, oo) 2 and any compact interval 
/ C (0, oo), we can obtain the convergence of generating functions of multitime correlation 
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functions, £^ n '°' L ' [x] — > G$[x]i as L -> oo. It implies p^ n ' 0,i ^ — y 3 p? as L — > oo in the sense 
of finite dimensional distributions weakly in the vague topology. In this case, we say that 
the noncolliding BESQ^- 1 (H^(t),P|) with an infinite number of particles £(R+) = oo is 
well defined with the correlation kernel [TO]. We will give sufficient conditions so that 
the process (Ey)(t), P£) is well defined, in which the correlation kernel is generally expressed 
using a double integral of an entire function represented by the Weierstrass canonical prod- 
uct having zeros on supp £, where supp £ = {x £ R : £({x}) > 0} (Theorem 12. 2p . As an 
application of this theorem, we will study the following example of infinite particle system, 
which is a non-equilibrium dynamics exhibiting a relaxation phenomenon. 
Consider the Bessel function 

J 'W = E r (B + ir( B + i + ,) (2) ' (116) 

n=0 

It is an analytic function of z in a cut plane. The function J u (z)/z u is an entire function. As 
usual we define z v to be exp(z/logz), where the argument of z is given its principal value; 



z = exp 



z/j log \z\ + \f— Targ(z)| , — 7r < arg(z) < it. (1-17) 
The function Jj,(z) is analytically continued outside this range of arg(z) so that the relation 

J v (e m7rV ^z) = e vm ^ l J v {z) (1.18) 

holds [29]. If v > —1, Jj,(z) has an infinite number of pairs of positive and negative zeros 
with the same absolute value, which are all simple. We write the positive zeros of J u (z) 
arranged in ascending order of the absolute value as 

< j v ,i < 3v,2 < jvfl <■■■ ■ 

Explicitly J u (z) is expressed using the infinite product of the Weierstrass primary factors of 
genus zero as (see Chapter XV of 



Jjz) 



[z/2f ff/ z 2 



The configuration in which every point of the square of positive zero of J u (z) is occupied by 
one particle, denoted by 

00 

i=i 

satisfies the conditions of Theorem 12.21 We will determine the correlation kernel of the 
noncolliding BESQ^ starting from explicitly (Theorem 12.31 (i)) and prove that the 
process shows a relaxation phenomenon to a stationary process, 

(SM(t + 0),P^)^(sW(t),PjJ as ^00 
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weakly in the sense of finite dimensional distributions (Theorem 12.31 (ii)). Here (S^(t), Pj v 
is the equilibrium dynamics, which is determinantal with the correlation kernel 



Kj„(t - s,y\x) = < 



due- 2u{s ~ t) J u {2^)J u (2 y /Tn]) if s<t 

J v (2y/x)jyJ'„(2jy) - ^J' v {2^)J v {2^y) 



x-y 

du e- 2u{s - t] J u {2^x)J u {2^y) 



if t 



if s > t, 



;i.2i) 



(x, y) E (0, oo) 2 , where J'(z) = dJ(z)/dz. Note that this kernel is temporally homogeneous 
but spatially inhomogeneous. Let be the determinantal (Fermion) point process on M. + , 
in which, for any iV E N, Appoint correlation function is given by [9] 



Pu{X N ) = 

with the Bessel kernel 

K Jv {y\x) 



det 

l<i,j<N L 



Kj v {xi\Xj / 



x N = (xi, . . . , xn) E (0, oo 



Kj v (0,v\x) 

J v (2y/x-)^yJ> v (2^y) - ^J' u {2^)J v {2^y) 



^J u+1 (2^)J u (2^y) - J v (2^x~)^yJ y+1 (2^y) 
x-y 

In particular, the density of particle at x E (0, oo) is given by 



;i.22) 



p v [x) 



K Jv {x\x) = lim K Jv (y\x) 

y->x 

(^(2v^)) 2 +(l-^) a(2v^)) 2 - 
(J,(2v^)) 2 - J v+1 {2y/x)J v . x {2y/x). 



;i.23) 



The probability measure is obtained in an A" — )■ oo limit called the hard- edge scaling- 
limit of the distribution of squares of eigenvalues of random matrices in the chiral Gaussian 
unitary ensemble [231 IS]- The process (E^(t),Pj u ) is a reversible process with respect to 
[A j v . The correlation kernel fll.2ip is called the extended Bessel kernel [TO], 127]. 

In the random matrix theory, three kinds of determinantal point processes of infinite 
particle systems have been well studied, the correlation kernels of which are given by (i) the 
sine kernel, (ii) the Airy kernel, and (iii) the Bessel kernel (231 IS]- They are obtained by 
taking (i) the bulk, (ii) the soft-edge, and (iii) the hard-edge scaling limits in the Gaussian 
unitary ensemble (GUE) for (i) and (ii) and in the chiral GUE for (iii), respectively. These 
three determinantal point processes have been extended to time-dependent versions so that 
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they describe equilibrium dynamics, which are reversible with respect to the determinantal 
point processes \T7\. Corresponding to these three stationary processes, the present authors 
introduced three relaxation processes with infinite numbers of particles realized in (i) the 
Dyson model (the noncolliding Brownian motion) starting from Z (i.e. the zeros of sin(7T2:)) 
in [19], (ii) the Dyson model with drift terms starting from the Airy zeros in [IB] , and (iii) 
the noncolliding BESQ^ starting from the squares of positive zeros of J„ in the present 
paper. The scaling limits are performed by increasing the number of particles in the system 
N — > oo, while in our setting determinantal processes with infinite numbers of particles 
N = oo are well constructed based on the theory of entire functions. In our relaxation 
processes in non-equilibrium we only have to wait for sufficiently long time to observe the 
three determinantal point processes. In order to give temporally inhomogeneous correlation 
kernels explicitly, we have reported the relaxation processes with the special initial config- 
urations. The universality of the three determinantal point processes in a wide variety of 
fields of mathematics, physics, and others [23l |9] implies robustness of relaxation phenomena 
with respect to initial configurations. Mathematical justification of this fact will be reported 
in the future. 

The present paper is organized as follows. In Section 2 preliminaries and main results are 
given. In Section 3 the properties of special functions used in this paper are given. Section 
4 is devoted to proofs of results. 

2 Preliminaries and Main Results 

We introduce the following operations; for £(•) = J2 i( z A 8 Xi (-) G 9Jt, 
(shift) with u E K, r^(-) = £/W(-)> 

ieA 



(dilatation) with c > 0, co£(-) = ^^S CXi (-), 
(square) £ <2) (-) = E^ 



■t6A 



ieA 



andfor£(-) = E ie AU-)enJT+, 

(square root) £<^>(.) = £ {s^ + S^)). 



iGA 



Note that the notation (11.201) states that this configuration is obtained as the square of 
the point-mass distribution on the positive zeros of J„(z) denoted by = X^i^ViO- 

We use the convention such that 



IJ/(x) = exp [ f Z(dx)logf(x)\ = J] /( 



x , €(W) 
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for £ G QJt and a function / on R. For a multivariate symmetric function g we write g((a;) x . g ,t) 
for g((xi) ieA ). 

The transition probability density of BESQ^, given by (II. 4p . satisfies the Chapman- 
Kolmogorov equation 

poo 

dypM(t-s,z\y)pM{s,y\x)=pM(t,z\x), (2.1) 

for < s < t, x, z G R+. Here we define the modified Bessel function of the first kind on C 

as 

r / \ / e— ^/ 2 J,(e^/ 2 z), -7T < arg(z) < tt/2, 

| ^^1/2^-3^/2^ < arg(2) < vr, 1 J 

where is defined by (I1.16P so that fll. 18j) holds [2J. This definition is consistent with (11. 5p . 
associated with the relation 

I u (e mn ^z) = e v ™^~ l l v {z). (2.3) 

For t G R, we define 



p M (t,?/|x) = <( 



1 (y\ v l 2 ( x + y\ ( \fxy s 

V — e~ y/2t , teR\{0},x = 0, 



(!) ^ exp (~^r) J * (^r) ' * G R \ {0} ' xe c \ {0} ' 



2|t| Vx, 

(2ltj)^r(i/ + i)' 

5(y — x), t = 0, x G C, 



(2.4) 

y £ C Then by using the modified version of Weber's integral of the Bessel functions 
(see Eqs.f l3.ip and (13. 2 j) in Section 3.1), (12. ip is extended to the following equations. For 

< s < t, x, z G C, 

/•oo 

dyp M (-t,z|?/)p M (t-s,?/|x) =p M (-s,z|x), (2.5) 
dyp {v \t- s,z\y)p {u) (-t,y\x) = p {u) (-s, z\x). (2.6) 

We define ^ 

( (-) P (,/ \t,y\x), teR,xeC\{0}, 
P jAt,y\x)=l W (2.7) 

[ y-^V^^^IO), t G R, x = 0, 

y G C. When t > 0, x, y G R + , it has the expression 



p Ju (t,y\x) = / duJ I/ (2v^)J„(2 v ^y)e- 2ut 
>/ o 

/■oo 

= 2 / dwwJ v {2w^x)J v {2w^/y)e- 2w2t . (2.8) 
9 



The Chapman- Kolmogorov equation (12.11) and its extensions (12. 5p and (12. 6p are mapped to 

dypj u (t - s,z\y)p Jv (s,y\x) = p Jv (t,z\x), (2.9) 

dypj v (-t,z\y)p Jv (t - s,y\x) =pj„(-s,z\x), (2.10) 
dy p Jv (t - s, z\y)p Ju (-t,y\x) =p Jv {-s,z\x) (2.11) 



for < s < t, x, z G C. 

For £ N e with fjv(]R + ) = iV G N, we define the functions of z G C, 

Ms*,*) = n ( 2 - 12 ) 

xeCivn{o} c 

n M (^,z) = ^ /2 n (£^). (2.13) 

For a G C we also define 

$o(6v,a,20 = n (r_ a e7v,^-a)= J] M - £^ J , (2.14) 



xeC]vn{a} c 



*M(6r,a,z) = <| W (2.15) 

nM(^,z), if o = 0. 

Theorem 2.1 (i) For any fixed configuration £ N G u»^^(M+) = iV G N, (SM(t),P£") 
zs determinantal with the correlation kernel 

Ki N (s,x;t,y) = lim— ?=/ dy' / tfepM(s, x|z)— ^$ &v, z, y')p {u) (-t, y'\y) 

-l(s>t)pM(s-t,x\y), (s,t) G [0,oo) 2 ,(x,y) G (0, oo) 2 , (2.16) 

where IV^jv) denotes a counterclockwise contour on the complex plane C encircling the 
points in supp £jy cn M + noi #ie point y' G (— oo, —e], and l(u>) is the indicator function 
of condition u>. 

(ii) If £n G 2^o" £iv(^+) =iVeN, £/ie correlation kernel is given by 

poo p0 

Kl»(s,x;t,y) = / £ N {dx') / dy'p^\s, x\x')%(£ N , x', y')p {u \-t, y'\y) 



-l{s>t)pM{s-t,x\y), (s,t) G [0,oo) 2 ,(x,y) G (0,oo) 2 . (2.17) 

Without changing any finite dimensional distributions of the process, the correlation kernel 
\2.17\) can be replaced by 

poo p0 

K j!,(s,x;t,y) = / £ N (dx') / dy' pj„(s, 2|x')$ H (6v, x', y')pj v (-t, y'\y) 

JO J-oo 

-l(s>t)p Ju (s-t,x\y), (s,t) G [0,oo) 2 ,(x,y) G (0, oo) 2 . (2.18) 
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Remark 1. If we consider a spatial distribution of particles at a single time t > of the 
noncolliding BESQ^, we have a determinantal point process with the correlation function 
Kl N (t,x;t, y), (x, y) E (0, oo) 2 . In particular, if we set s = t = 1/2, the correlation kernel 
(I2.16P is reduced to be the kernel of the perturbed chiral GUE of random matrices given in 
Proposition 5 by Desrosiers and Forrester [8]. More detail, see Remark 2 in Section 3.3. 

For L > 0, a > and £ G 9JI we put 



M a {t,L) 
and 



■L,L}\{0} \ x \ a 



M a {i)= lim M a (£,L), 

L— >oo 

if the limit finitely exists. We introduce the following conditions for configurations £ G 50T + . 

(C.A) (i) There exists a E (1/2, 1) and C x > such that Af a (f) < Ci. 
(ii) There exist /3 > and C2 > such that 

M^r^O < C 2 (|o| V 1)-?, VaG supp£. 



We denote by X + the set of configurations satisfying the conditions (C.A), and put Xq = 
X + fl 9Jto . For £ G Xq, a G ffi. and z G C we can define 

$o(£, a, z) — hm ^o(£ H [a — L, a + L], a, z), 

L— >oo 

a, z) = lim $ M (£n [a - L, a + L], a, 2), 1/ > -1. 

L— >oo 

Since $o(£ ; 2) has the expression of the Weierstrass canonical product of genus zero (see 
(12.141) ). it is an entire function of a variable z G C [22]. The set of zeros of a, 2) is 

given by (supp £U{0})fl{a} c and all zeros except are simple for £ G 971q . 

Theorem 2.2 If £ E the process (H^(t), P£) zs we// defined with the correlation kernel 

poo r0 

Ki(s,x;t,y) = / £(aV) / dy' p w (s, x|x')$o(6, 2/> H (-*, !/'|y) 

JO J -oo 

-l(8>t)pM(8-t,x\y), (s,t) E [0,oc) 2 , (a;,?/) G (0,oo) 2 . (2.19) 
This correlation kernel 1K| can be replaced by 

/■oo /"0 

K* Jv {s,x;t,y) = / £(aV) / dy'pj„(s, x\x')&»\t, x', y')p Jv {-t, y'\y) 

J0 J -00 

-l{s>t)p Jv (s-t,x\y), (s, t) E [0, oo) 2 , (x, y) E (0, oo) 2 (2.20) 
without changing any finite dimensional distributions of the process. 
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In case £(K+) = oo, Theorem 12.21 gives the noncolliding BESQ^ with an infinite number of 
particles starting from the configuration £ £ Xq. It is easy to check that, if (x, y) £ (0, oo) 2 , 

Kl(t,x;t, y)Kl(t,y;t, x)dxdy — >■ £(dx)l(x — y), t — V in the vague topology. 

For 7 > 0, we put 

g 1 (x) = x 7 , x £ R + , 

and 

oo 

^(•) = Ew-)- 
i=i 

For any 7 > 1 we can show by simple calculation that rf 1 satisfies (C.A)(i) with any 
a £ (I/7, 1) and some C\ = C\(a) > depending on a, and (C.A) (ii) with any f3 £ (0, 7 — 1) 
and some Ci = C*2(/3) > depending on (3. This implies that rp is an element of Xq for any 
7 > 1. 

More interesting example is given by the following theorem. 

Theorem 2.3 (i) The noncolliding BESQ" starting from Qj , (S M (t), K Jv ), is well defined 
with the correlation kernel 

^J„(s,x;t,y) = jr f dzpj^s^Uli)-^ 1 — f" * . pj v (-t,z\y) 

-l(s > t)p Jv (s - t, x\y), (s, t) £ [0, oo) 2 , (x, y) £ (0, oo) 2 . (2.21) 

(ii) Let (E^(t),P j v ) be the equilibrium dynamics, which is determinantal with the extended 
Bessel kernel M.21\) . Then, for t > 

(E^(t + 9),F?)^(E^(t),PjJ as 9^oo (2.22) 
weakly in the sense of finite dimensional distributions. 

3 Some Properties of Special Functions 

3.1 Integral formulas of Bessel functions 

The following integral formulas are known [291 E]- For $ti> > —l,p, a, b > 0, 

J~ 'duue-^ 2 J u (au)Ubu) = ^exp (~^) h (0) , (3.1) 
[ duue-^Uau)Ubu) = ± exp I v (0) . (3.2) 
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The equalities flU}, ((23]) and ([2l| for and f[2T8]) - f[2~TTj) for p Ju are derived from the 
above integral formulas. In addition to them, the following equality is also derived from 
f!3~Tj) with (P22D and (|2Tg]l . For t > 0, x, z > 



dy J u (y/zy)p Ju (-t,y\x) = e tz/2 J u {yfzx). 
Lemma 3.1 For any i E N,z ^ 

1 duuJ v {zu)J y {j v ^u) 



(3.3) 



jld ~ Z 2 Jv+l(jv,i) (Ju+x{jf,i)) 2 JO 



1 1 



(^+l(j^)) 2 2 7o 

Proof. The following formula is found on page 482 in [2 



dw J y (Zy/w) J u (j Uii y/w) . 



(3.4) 
(3.5) 



f 

■Jo 



tJ v (at)J v (aoi)dt 



x 



a' - «n L 



J u (ax)aoJl(aox) — J u (aox)aJl(ax) 



Set a = j V) i and x — 1. Since J u (jv,i) = by definition of j u /s, we have 



tJ v (at)J v (j Uti t)dt 



« 2 - J*J< 



Note that the Bessel function satisfies the relation J' v {z) = (v/z)J v {z) — J u+ i(z). Then 

■^OW) = - J v+iUv,i) ( 3 - 6 ) 

and thus 



tJ v (at)J v (J V)i t)dt 



? 2 • — a 2 



Divide the both sides by (J u +i(ju,i)) 2 , we obtain (13. 4p . Eq. f)3.5p is obtained from (13.4)) by 
setting k = y^. | 

3.2 Fourier-Bessel expansion 

As shown in Section 18.24 in any continuous function f(x) e £ 2 (0, 1) has the expansion 

2 



(4+l(j^)) 2 7o 



uf{u)J v {j V)i u)du. 



That is, 



/(*) 



f duuf(u)2_^ — — — — — — , feL(Q, 1) 



v+iOW)) 5 
13 



which is called the Fourier-Bessel expansion. Set u = ^Jy and then replace x by t/x, we 
have 

Ju(ju,iVx)Jv(jt>,iy/v) 



/(VE)= / dy/(Vy)X) 

^0 „-_1 



In other words, the functions { J v (j v ,i\/x) / Ju+i(ju,i), i G N} form an orthonormal basis for 
/ G L 2 (0, 1) and the completeness is also established; 

= %-!/). »,»6(0,1). (3.7) 

3.3 Multiple orthogonal polynomials 

Fix a configuration 

6v e Wl + with &v( K +) =iVeN. 

In the present paper the multiple orthogonal polynomials associated with the modified Bessel 
function l v indexed by £n are defined by the following 0, EH [8] . 

Type I: The multiple orthogonal polynomials of the type I are the set of functions 

^A^(y,x) : x G supp £jv, polynomial of y of degree £at(x) — 1 j (3.8) 

such that, if we set 

<4>) = £ < } (l/^)^(f)' /2 /,(v^)e- ( ^ )/2 , (3.9) 

xGsupp £iv 

then 

f *«S3toV={^ iZl^if^- 2 (3.io) 

Type II: The multiple orthogonal polynomial of the type II is the monic polynomial of degree 
6v(K+), 

p£ ) (y) = y Cw( ^ ) + 0(^>- 1 ) (3.11) 

such that for each x G supp £jv 

dyPtJ(y)y'-(l) h(V^y)^ ix+y)/2 = o, o<i<Ui)-i- (3.12) 

The following integral representations have been obtained by Desrosiers and Forrester 

0- 
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Lemma 3.2 The functions Q^(y) and P^(y) have the following integral representations , 



PtJ(y)= f dwU-Y I^-^e^'Hiw-x), (3.14) 



where r(£/v) denotes a counterclockwise contour on the complex plane C encircling the points 
in supp £tv on R + . 

Remark 2. The present definition of multiple orthogonal polynomials associated with the 
modified Bessel function I u is slightly different from that given by Desrosiers and Forrester 
[8]. Moreover, since they have used the function Fi(a;z) = Xl^Lo zU / {( a )n n ^} in order to 
express the polynomials instead of I v (z), our expressions (I3.13P and (I3.14p seem to be quite 
different from their functions. The identity 

is established, however, and then we can see that Lemma 13.21 given above is equivalent with 
Proposition 6 of [8J. More precisely speaking, if we write the orthogonal polynomials in [8] 
as Qti F (y) and P^ F (y), where the parameters a = u, n = (ni, . . . , no), and a = b n , we have 
the correspondence 

gg(y) = (I) , pM(y) = 2l«>? F (|) , (3.15) 

for = YaLi = Ye=i n ^2b e (-) and &v( R +) = N = J2?=i n £ = 

We write £/v(") = Si=i with a; = (a^)^ such that < x\ < %2 < • • ■ < xn- Then 
we set 

i 

eiv,o(-) = and eJV,<(0 = X)M0, l<^<iV- 

fc=l 

By this definition 6v,i(R+) = i,0 < i < N and ^({x}) < 6v,i+i({^}), Vx G R+,0 < z < 
N — 1. We define 

^ ,+) (y;eiv) = Qg i+1 (y), m^;£*) = p&>), o<<<jv-i. (s.ie) 

By the orthogonality relations (I3.10p . 03.121) and the above definitions, we can prove the 
biorthonormality 

POO 

/ dyMf>- ) (y;£ N )Ml v > + \y;Z N )=6 ihl < i, k < N - 1. (3.17) 
Jo 
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Lemma 3.3 Let N e N, £ N e M + with £ N (I 
N-l, 



dxp^(t-s,y\x)M^' +) l-;-o£ N 



s s 



N. For < s < t, x, y G R+, < i < 

S V +1 M (^+) fi. 1 



dyM\ 



(",-) 



X 1 

1 

s s 



°6v 

°6v 



(3.18) 
(3.19) 



Proof. By definition (13. 16p and Lemma [3.21 



(",+) 



x 1 



s s 



IN 



r((i/s)o^, i+1 ) 



, 1 1 /X\»P 

U>Z 777 — 

2z v ' 2 \s 



a/s) 



By setting z = w/s and by the definition (12.41) of p^ v \ we find 



M 



(",+) 



x 1 



s s 



2ttV 3 T 
1 



dw ■ 



1 /xW 2 



r (&v,t+i) 



2s Vw; 



^ e -( w+x )/(2s) 



n 



u> — a) 



r(£jv,i+i) 



27TV-1 

Therefore the LHS of (ETTBjl is 

dw 



dwp^ u \s, x\w) 



(3.20) 



J+i 



1 



T(^ ii+1 ) n a eCiv, i+1 ( w _ a ) Jo 
By the Chapman-Kolmogorov equation (I2.ip . it equals to 



dxp^(t — s, ?/|x)p^(s, x|w). 



1 



2^^ Jr(^, i+1 ) IW >i+1 (™ - a) 



Then, comparing with the expression (I3.20p . we have (I3.18p . 
By definition (I3.16P and Lemma l3~2l 

u/2 



By setting w = u/t and by the definition (12. 4p of , we have the expression 



M. 



dwp (l/) (-t,w|?/) J| (it -a). 



(3.21) 
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Therefore the LHS of fl3TT9j) is 

du II ("-<') / dyp (u \-t,u\y)p {iy \t- s,y\x). 

-oo „ c , JO 

By the extension of the Chapman-Kolmogorov equation (I2.5p . it equals to 

f-0 



f 



du Y[ (u-a)pM{-s,u\x)= (^Ym^'- ] 1 ' 



s s 



N 



where the expression (13.211) was used. This completes the proof of (I3.19p . | 

Lemma 3.4 Let y = {yi)f =1 G W^. For any £jv(-) = X)i=i e ^ + w ^ a labeled 

configuration x = (xi)f =l such that %\ < x<i < ■ ■ ■ < xn, 



, , \ det 

h N (x) l<i,j<N 



U|)' /2 e- ( " +%)/ %(v^: 



det 

l<i,j<N 



(3.22) 



Here, when some of the Xi's coincide, we interpret the LHS using VHopitaVs rule. 
Proof. By the multilinearity of determinant 



7 — 7 — r det 

h N {X) l<i,j<N 



1 ( Vj 



det 

l<i,j<N 



det 

l<i,j<N 



1 ( Vi 



2 \xi 



2 \ Xi 
u/2 



u/2 



-fo+tt)/2j f/^y- ) . 1 

llfc=lv J ' 1 



1 / 



En* 1 • 



i//2 



1 



= 1 - \— <-/ i il<k<i,kJ4\ x t Xk 

By definition flBTTHj) with f l3TT3|) of Lemma 13721 when &v G 97tf, &v(K+) = N, 



z — X) 



2 ^ 'nut*-** 



2W-1 



^2 la* 



u/2 



- {xe+y)/2 iu(V^y) 



IL 



<k<i,k^l 



(x e - X k ) 



(3.23) 



1 < i < N, y G R+. Then (I3.22|) is proved for £jy G SPTq". When some of the x^s coincide, the 



LHS of fl3T22|) is interpreted using l'Hopital's rule and in the RHS of (ET22jl 
should be given by the first expression of (13.231) . Then (I3.22j) is valid for any £ N G 97t + with 
6v(R + ) = JV e N. i 
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4 Proofs of Theorems 



4.1 Proof of Theorem 12.11 



In this subsection we give a proof of Theorem 12. 1L Assume that £n G with £tv(R + ) 
N EN. Define 



TV 



(4.1) 
(4.2) 



< i < N — 1, t > 0, x G R. From Lemma [3 .3[ the following relations are derived. 
Lemma 4.1 For £ N G 9H wift &v(R) = N G N, < t x < t 2 , 

dx 2 4"'~\t2, x 2 ; £ N )p iv) (h - t u x 2 \x 1 ) = 4 v '~\h, xi; 6v), < % < N - 1, (4.3) 
dx 1 pM(t 2 -t 1 ,x 2 \x 1 )<f>!i v ' + \t 1 ,x 1 iZ N ) = 4 u ' + \t 2 ,x 2 ^ N ), 0<i<N-l, (4.4) 

POO 

dxi / dx 2 (j)^'~\t 2 ,x 2 ;^ N )p^\t 2 -t 1 ,x 2 \xi)(j)f' + \t 1 ,Xi;^ N ) = S { 



< i, j < N - I. 



(4.5) 



Remark 3. The equation (14. 5 p is obtained from the combination of (14. 3 p and (14.41) . It 
should be emphasize the fact that the biorthonormality (13.171) is obtained by just taking the 
limit t 2 — t\ in (14.51) . For the purpose in this paper, the connection to the theory of 
multiple orthogonal polynomials [HI 0, El Ell E] is not so essential. We can think that the 
integrals (I3T2TJ and (EE2I1) define the functions (j)f & \t,x;£) and </>[ v '~\t,x;g) through (Q]l 
and (14.21) . respectively. As shown in the proof of Lemma |3.3[ (14. 3p and (14.41) are readily 
obtained by applying the Chapman-Kolmogorov equation (12. ip and its extension (12. 5p . 

By definitions (13.111) . (13.161) . and (14.21) . we can see that (jxf'~'(t, x; £n) is a monic polyno- 
mial of x of degree i, which is independent of t. By this fact and Lemma 13 .41 the multitime 
probability density (11.101) is expressed as 



r£ N (t, -t*r £ {M) ) - Het 

Vv \ l l> SAT j ■ ■ ■ i l M, VZV ) — , 



<Pi-l [tM,Xj ,C,N) 



M-1 



X 



n f [ N^ 



m+1 ^rai 



t m :x (m+1) \ X ( m) ) det 



m=l 



l<k,£<N 



0fct > (*i,4 1) ;6v) 



(4.6) 



for £ N G £0T + with £at(R + ) = N G N. By the argument given in Section 4 in [T7J, the 
expression (14.61) with Lemma 14.11 leads to the Fredholm determinantal expression for the 
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generating function of multitime correlation functions (I1.12p . 



< 31 N [X} = Det 

\<m,n<M 
(a;,i/)e(0,oo) 2 



S m J{x -y) + S m ' n (x, y; £ N )Xn{y) 



where 
with 



S m ' n (x,y;£ N ) = S m > n (x,y;£ N ) - l(m > n)p^(t m - t n ,x\y) 



N-l 



L m \ b m / \<"m »m / \ b n b n 

Here the Fredholm determinant is defined by the following expansion 



o£ N . (4.7) 



Det 

l<m,n<M 
(z,y)e(0,oo) 2 

2V N M 



S mn S{x -y) + S m ' n (x, y; £ N )xn{y) 



N N M JVi r N M 

Vi=0 AT M =0m=l m " /K + i=l i=l 



iVi=0 A^Af 

M AT, 

X 



nri*. 



m \ Xj J det 

l<i<iV TO ,l<i<^Vn 
m=1 4=1 l<m,n<M 



(4- 



The following invariance of finite dimensional distributions of determinantal processes 
will be used. 

Lemma 4.2 Let (H(t),F*) and (H(t),P ? ) be the m + -valued processes, which are determi- 
nantal with correlation kernels K and IK, respectively. If there is a function G(s,x), which 
is continuous with respect to x G (0, oo) for any fixed s G [0, oo), such that 



K(s,x;t,y) 



G(s, x) 



K(s, x; t, y), (s, t) G [0, oo) 2 , (x, y) G (0, oo) 1 



then 



(H(f),P*) = (3(f), 
zn i/ie sense of finite dimensional distributions. 



(4.9) 
(4.10) 



The relation (14. 9p is called the gauge transformation and (14. 10p is said to be the gauge 
invariance of the determinantal processes 
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Proof of Theorem \2.1[ Inserting the integral representations for M- u,± \ given by ( I3.20p and 
(HL2D, into dUTJ), the kernel S m > n is written as 



1 '° 



S m > n (x,y;£ N ) = — = / du & dzpM(t m ,x\z)pM(-t n ,u\y) 

n?=i( M - Q j) 
fc=0 n i= iu-%) 

where ^jv(-) = ^oi(')- For ^Ij 2 ^ £ C with z\ ^ {xi, . . . , xn}, the following identity 

holds, 

By this identity, we have 

27TV-1 J-oo JVuit^o^) 



N 

U — Z \ -fJj 2 — Oj 



Note that 

1 /" _e /" 1 

du (p dzp {u \t m) x\z)p {u \-t n ,u\y) = 0, \/e > 0, 

it — z 



27TV-1 7-00 h u {f^: i N ) 

since by definition (12.41) 



P {t m} x\z) {2tm)v+x e 2^T(n + l)T(n + l + u)(2t mJ 



\2n 



is an entire function with respect to z. Then (I2.16P is obtained. When £jv G , the Cauchy 
integrals are performed in (12 . 1 6[) and (12.171) is obtained. If we use (\2.7\i and (12.151) . we have 
the equality 

fx \ v/2 

Ki»(s,x;t,y) = - Kf%(s,x;t,y), (s,t) e [0, oo) 2 , (x, y) e (0,oo) 2 . (4.11) 

Then, by Lemma l4T2| the correlation kernel can be replaced by IKj^ without changing 
any finite dimensional distributions. This completes the proof. | 
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4.2 Proof of Theorem 12^21 

In this subsection we give a proof of Theorem 12.21 For £ G QJt, L > 0, we define 



M(C,L) 



[-l,l]\{o} a; 



and put 



M(C)= hm M(£,L) 



if the limit finitely exists. In an earlier paper we proved the following lemma (Lemma 4.4 in 

USD- 

Lemma 4.3 Assume that £ G QJt satisfies the following conditions: 
(C.l) t/iere exzsis C > smc/i i/iai |M(C)| < Co, 

(C.2) (i) t/iere exzsi a G (1, 2) and C\ > snc/i #ia£ Ms(£) < Ci, (ii) t/iere exzsi /3 > and 
C2 > suc/i that 

M!(r_ a 2C <2> ) < a 2 (|a| V Va G supp C 

TTien t/iere erases C3 = Cs(a, (3, Co, Ci, C2) > and 6 G (a V (2 — 2) smc/i i/ia£ 

, G K, Va G supp £■ 



|$o(C, a, < [C 3 {(M V If + (|a| V if 

Here we use the following version. 

Lemma 4.4 For any £ G j£q , there exist C3 = Cs(a, (3, Ci, C2) > and 9 G (a V (1 — /3), 1) 
sitc/i t/iat 

|d>o(£,a,y)|<exp [C 3 {(|a| Vl) 9 + (|y| Vi; 



, Vy < 0, Va G supp £. 
Proof. First we consider the case that a = G £. Let £ = £( 1//2 >_ Then the equalities 



n i' 

xG?n{o} 

n 11 

x'GCn{o} 



n " 

xG^n{o} 
= $o(C,0,v^) 



1 + 



(4.12) 



hold, if the products finitely exist. Since M(Q = by the definition £ = £^ 1//2 \ C satisfies 
the condition (C.l) in Lemma [4. 31 By assumption of the present lemma, that is, £ satisfies 
the condition (C.A), the condition (C.2) is satisfied by £ with a = 2a and (3 = 2(3, since 

M a (0 = (M 2a (C)) 2 and M x (r_ a = M x {t_ { ^ 2) ), a G supp £. 

Then, by Lemma S31 there exist G (S V (2 - 0), 2) = (2a V 2(1 - (3), 2) and C 3 < 00 such 
that <3>o(C 0, v 7 -^) < exp C^^fz VI) 6 * , z > 0. Since ( I4.12p is established, it implies 



$o(e,0,-z)<exp C 3 (|z|Vl) 



(4.13) 
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with 9 = 9/2. 

Next we consider the case that a G supp £ and a / 0. For y < 0, the equality 
$ (£, a, y) = $ (e, 0, y)$ (£ H {0} c , a, 0) (^J 



a- y 



is valid. Since y < is assumed 



y\?«o}) a 
aJ 



On the other hand, 

$o(£n{0} c ,a,0) 



a - y 



n 

ae£n{o} c 



a-y 



V 



2/ 



a - y 



< 1. 



1 + 



x — a 



< 



expjy (r_ £)(da;)|^|| = exp ||a|Mi(r_ a £)|. 



By the condition (C.A) (ii) of jEq , it is bounded from the above by exp{C 2 (|a| V 1) P\a\} = 
exp{C*2(|a| V 1) 1_/3 }. Combining this with ( I4.13p . the proof is completed. | 

Proof of Theorem {EM Note that £ n [0, L], L > and £ satisfy (C.A) with the same 
constants C^C^ and indices a, (3. By Lemma [4.41 we see that there exists C3 > such that 



|$ (en [0,L],x',y')| < exp C 3 (|x'| V l) 9 + (|y'| e V 1) 



VX > 0, W G supp £ n [0,L], Vy' < with 9 E (a V (1 — 1). Therefore, since for any 
x' G supp £, y' < 

$ (£ n [0, L], x', y') -> $o(C, x', 2/'), L -> 00, 
we can apply Lebesgue's convergence theorem to (12.181) and obtain 

lim Kf 1 ^ ( S ,x;t,y) = K« ( S ,x;t,y). 
Since for any (s,t) G (0, oo) 2 and any compact interval I C (0, 00) 



sup 



< OO, 



we can obtain the convergence of generating functions for multitime correlation functions, 
<?^°' L '[xl ~~ y @tix], as L — )• 00. It implies p-^ '^ — > as L — > 00 in the sense of finite 
dimensional distributions. By virtue of the gauge invariance of determinantal processes 
(Lemma 14.21) . the proof is completed. 1 
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4.3 Proof of Theorem Q 

In this subsection we give a proof of Theorem 12.31 First we prove a lemma. 
Lemma 4.5 For z ^ j Ut i, z G N 

*"^™-u$-*'Xr (4 ' 14) 

Proof. By the expression flEDg) of J u (z) and the definition (I27T5ji of B», we have 

= f^hf M{ ^ y) - (415) 

Put nM'(e <2> ,^ 2 ) = (d/dz)U^(^ 2 \z 2 ). Then 

= f^T) nH/ (^' &)• ( 4 - 16 ) 



We see that, if x' G £, 



xe£n{z'} c 

Since 



if a;' G f , 

nH' ( ^,(xo 2 )=-5(xT n f 1 - 1 ^) 



Therefore, for x' G £, 

i n^V 2 ) _*V£\" n i - ^ 2 /^ 2 r417 s 

l-z 2 /(x') 2 (UM)'(^ 2 ),(x') 2 ) 2\x'J J-J- l-(x') 2 /x 2 ' l ' ' 



xe£n{x'} 

Since 



z 2 — a 2 x 2 — z 2 1 — ^ 2 /x 2 



x 2 — a 2 x 2 — a 2 l — a 2 /x 2 
for x a, combining (14. 17p with the definition ( I2.15j) gives 

Z\ v t-i- 1 — 2 2 /x 2 



$H(e^,(xo 2 ,, 2 ) = (^y n i 



,x') 2 /x 2 

x' G £. (4.18) 



^ _ ( X ') 2 nH'(e< 2 ),(x') 2 ) 
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Setting £ = and x' = j U;i in (j4.18j) . we have 



by (OHjl and (OB]l . If we use (M . (Oljl is obtained. i 

Proof of Theorem \2.<A We set £ = £jf in and obtain the kernel 

oo „o 

K^( S ,x;t,y) = dzp Jv (s,x\jl t ^ v \^,Uu^ 2 ,z) P j u {-t,z\y) 
i=i 

-l(s>t)p Jw (s-t,ar|3/). (4.19) 

By Lemmas 13.11 and I4.5[ 

J v {\fz) 



Uu,i) 2 - Z J u +l{ju,i) 

du J v ( y/uz) Jy ( ) • 



By the integral formula fl3.3|) . the first term of the RHS of (14.191) equals 

y^J„(s,£|(j^) 2 )— — u 2 / due ut/2 J v (yfmj)J v (y/uj^i). 
Since s > 0, we can use the expression (12. 8 p for Pj v (s, x\j 2 J and the above is written as 

Jo Jo ^ (^+i(>,i)) 2 



From the completeness (13.71) of the orthonormal basis {J v (jv,i\/x) / 'J v +i(j V j),i G N} in 
L 2 (0, 1), the above gives 

Kj u (s,x;t,y) = K Ju (s,x;t,y) + R(s,x;t,y) 
with the extended Bessel kernel (ll.2ip and 

pl POO 

R(s,x;t,y) = du dw e ut/2 - 2ws J u (y/uy) J u (2^Jx) 



o Jl 

oo 



E J V (2y/wj l/i i)J u (y/u j^ i ) 
( T Jn 



Since for any fixed s, t > 

|i2(s + 0,a; : t + 0,y)| as 6> -» oo 
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uniformly on any compact subset of (x,y) <G R+, 

K Jv (s + 9,x;t + 9) ->K Jv (s,x;t,y) as 9 -> oo 

holds in the same sense. Hence we obtain fj2. 22(1 . This completes the proof. | 
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